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ABSTRACT 

Progenitor scenarios for short gamma-ray bursts (short GRBs) include coalescenses of two neutron stars or a 
neutron star and black hole, which would necessarily be accompanied by the emission of strong gravitational 
waves. We present a search for these known gravitational- wave signatures in temporal and directional coincidence 
with 22 GRBs that had sufficient gravitational-wave data available in multiple instruments during LIGO's fifth 
science run, S5, and Virgo's first science run, VSR1. We find no statistically significant gravitational-wave 
candidates within a [—5, +1) s window around the trigger time of any GRB. Using the Wilcoxon-Mann- 
Whitney U test, we find no evidence for an excess of weak gravitational-wave signals in our sample of GRBs. 
We exclude neutron star-black hole progenitors to a median 90% confidence exclusion distance of 6.7 Mpc. 
Subject headings: compact object mergers - gamma-ray bursts - gravitational waves 
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1. INTRODUCTION 

The past decade has seen dramatic progress in the under- 
standing of gamma-ray bursts (GRBs), intense flashes of 7- 
rays that are ob served to be isotropically distribute d over the 
sky (see e.g. Klebesadel et al. 1973 Meszaros 2006 and refer- 
ences therein). The short-time variability of the bursts indicates 
that the sources are very compact. They are observed directly 
by 7-ray and X-ray satellites in the Interplanetary Network 
(IPN) such as HETE, Swift, Konus-Wind, INTEGRAL, and 



primarily in the coalescenc e of a neutr on star (NS) with another 
compact object (see, e.g., Nakar 2007 and references therein), 
such as a neutron star or black hole (BH). There is growing 
evidence that finer distinctions may be d rawn between bursts 



Ferm i (see|Ricker et al.[2 003 ; Geh rels et al.[ 2004 ; Aptekar et al. 
[l9l^|Wmkler et al.|2003HAtwood et al.|2009| and references 

therein). 

GRBs are usually divided in to two types (see |Kou veliotou | 
|et al.|1 993 ; Geh rels et al.|20 06 ), distinguished primarily by the 
duration of the prompt burst. Long-duration bursts, with a du- 
ration of > 2 s, are generally associated with hypernova explo- 
sions in star-forming galaxies. Several nearby long GRBs have 
been spatially and temporally coincide nt with core-collapse 
supernovae as observed in the optical ( Campana et al.|2006| 



Gala maet"aT.| 1998) |Hjorth et al.|2003| pvlalesani et al.|2004| ). 
Follow-up observations by X-ray, optical, and radio telescopes 
of the sky near GRBs have yielded detailed measurements 
of afterglows from more than 500 GRBs to date; some of 
these observations resulted in strong host galaxy candidates, 
which allowed redshift determination for more than 200 bursts 
( |Greiner|2009l >. 
Short GRBs, with a duration < 2 s, are thought to originate 
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(Zhang et al. 2007 Bloom et al.||2008[ ); for example, it is 
estimated that up to ~15% of ihort GRBs could b e associated 
with s oft gamma repeaters ( Nakar et al.|2006 Chapman et al. 
2009), which emit bursts of X-rays and gamma rays at irregular 



intervals with lower fluence than co mpact binary coalescence 
engines ( |Hurley et al.|2005||pilrner et al.|2005) . 

In the compact binary coalescence model of short GRBs, a 
neutron star and compact companion in otherwise stable orbit 
lose energy to gravitational waves and inspiral. The neutron 
star(s) tidally disrupt shortly before coalescence, providing 
matter, some of which is ejected in relativistic jets. The prompt 
7-ray emission is widely thought to be created by internal 
shocks, the interaction of outgoing matter shells at different 
velocities, while the afterglow is thought to be created by 
external shocks, the inter action of the outflowing m atter with 
the interstellar medium ( Meszaros 2006) |Nakar|2007 1 ) . If the 
speed of gravitational radiation equals the speed of light as 
we expect, then for an observer in the cone of the collimated 
outflow, the gravitational-wave inspiral signal will arrive a few 
seconds before the electromagnetic signal from internal shocks. 
Several semi-analytical calculations of the final stages of a NS- 
BH inspiral sho w that the majority of matter plunges onto the 
BH within 1 s ([Dayjes et al. 2005). Numerical simulati ons on 
the mass transfer su ggest a timescales of milliseco nds (|Shibata| 
|& Taniguchi]2 008) or some seconds at maximum (Faber et al. 
2006b. Also, it has been found in simulations that the vast 
majority of the NS matter is accreted onto the BH directly and 
promptly (with in hundreds of ms) without a toru s that gets 
accreted later (Rosswog 2006^ |Etienne et al.|2008| l. 

Compact binary coalescence is anticipated to generate strong 
gravitational waves in the sensitive frequency ba nd of Earth- 
based gravitational-wave detectors (Thorne 1987). The direct 
detection of gravitational waves associated with a short GRB 
would provide direct evidence that the progenitor is indeed 
a compact binary; with such a det ection it would be possi- 
ble to measure compon ent masses ( [Cutler & Flanagan|1994[ 
Finn & Chernoff 1993), measure component sp ins (fPo isson 
& Will||1995|l, constrain NS equa tions of state ( {Flanagan & 
Hinderer||2008[ |Read et al.| 2009), test general relativity in 
the strong-field regime (Will 2005), and measure calibration- 
free luminosity distance (Nissank e et al.|[20 09), which is a 
measurement of the Hubble expansion and dark energy. 

In this paper, we report on a search for gravitational-wave 
inspiral signals associated with the short GRBs that occurred 
during the fifth science run (S5) of LIGO JAbbott et alT|2 009a), 
from 4 November 2005 to 30 September 2007, an d the first 
science run (VSR1) of Virgo ( jAcernesee t al. 2008 ), from 18 
May 2007 to 30 September 2007. S5 represents the combined 
operation of the three LIGO detectors, one Michelson inter- 
ferometer with 4 km long orthogonal arms at Livingston, LA, 
USA, named LI, and two interferometers located at Hanford, 
WA, USA, named HI and H2, with lengths of 4 km and 2 km 
respectively. VSR1 represents the operation of the Virgo in- 
terferometer located at Cascina, Italy, named V 1 , which has a 
length of 3 km. During the S5/VSR1 joint run, 212 GRBs were 
discovered by different satellite missions (39 of them during 
VSR1 times), 33 of which we clas sified as search targets (8 of 
them in VSR1 times). See section [2/2] for more details on the 
selected GRBs. 

A similar search in the same LIGO/Virgo data-set was per- 
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formed in Abbott et al. ( 2009c| l, looking for short-duration 
gravitational-wave bursts in association with 137 GRBs 
recorded during S5/VSR1, both long and short. The anal- 
ysis reported upper limits on the strain of a generic burst of 
circularly polarized gravitational radiation, predominantly at 
the detectors' most sensitive frequencies. These were trans- 
lated into lower limits in distance by assuming that 0.01 M 
is converted into isotropically emitted gravitational waves. In 
contrast, the search described in this paper does not make any 
assumption on the polarization of the gravitational waves and 
searches for the specific signals expected from binary coa- 
lescenses. Importantly, the present search can distinguish a 
coalescence signal from other models and estimate the progen- 
itor parameters. 

The remainder of the paper is organized as follows. In Sec- 
tion^ we discuss the set of GRBs we chose for this analysis 
and outline our analysis methods. In Section [3] we present 
the results and astrophysical implications for the GRBs in our 
sample. 
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2. SEARCH METHODS 

2.1. Experimental setup 

The binary coalescence model suggests that the time delay 
between the arrival of a gravitational wave and the arrival of 
the subsequent electromagnetic burst, referred to as trigger 
time, is a few seconds. We assessed uncertainties in reported 
trigger times and quantization in our own analysis along integer 
second boundaries, finding that these each contribute less than 
one second. For example, when the Swift BAT instrument 
determines that the count rate has risen above a threshold, 
it waits for the maximum to pass, checking with a 320 ms 
cadence (Gehrels et al. 2008); it reports the start time of the 
block containing the maximum, rather than making any attempt 
to identify the start of the burst, and does so with a 320 ms 
granularity. As another example, there h ave been reports of 
sub-threshold precursors to many GRBs (Burlon et al. 2009 



For each GRB in our sample, we checked tens of seconds 
of light-curve by eye to look for both excessive difference 
between the trigger time and the apparent rise time, and also for 
precursors, but found nothing to suggest that we should correct 
the published trigger times. The largest timing uncertainty 
we identified is the delay between the compact merger and 
the prompt emission of the internal shocks. We search for 
gravitational-wave signals within an on-source segment of 
[—5, +1) s around each trigger time for each GRB of interest, 
feeling that this window captures the physical model with some 
tolerance for its uncertainties. 

Because we believe that a gravitational wave associated with 
a GRB only occurs in the on-source segment, we use 324 
off-source trials, each 6 s long, to estimate the distribution of 
background due to the accidental coincidences of noise triggers. 
We also re-analyze the off-source trials with simulated signals 
added to the data to test the response of our search to signals; 
these we call injection trials. The actual number of off-source 
trials included in the analysis varied by GRB, as the tri als that 
overlapped with data-quality vetoes were discarded (Abbott 
|et al.|2009d[ ). To prevent biasing our background estimation 
due to a potential loud signal in the on-source trial, the off- 
source segments do not use data within 48 s of the on-source 
segment, reflecting the longest duration of templates in our 
bank. Finally, we discard 72 s of data subject to filter transients 
on both ends of the off-source region. Taking all of these 
requirements into account, the minimum analyzable time is 



Table 1 

Parameters of the 22 GRBs selected for this search. The values in the 
References column give the number of the GRB Coordin ates Network (GC N) 
notice from which we took the preceding information iBarthelmy 20091. 

2190 s. 

2.2. Sample selection 

X-ray and 7-ray instruments identified a total of 212 GRBs 
during the S5 run, 211 with measured durations; 30 of them 
have a T90 duration smaller than 2 seconds. T90 is the time in- 
terval over which 90% of all counts from a GRB are recorded. 
While the T 90 classifies a burst as long or short, it is not a 
definitive descriminator of progeni tor systems. In ad dition 
to the short G RBs, GRB 051211 ( |Kawai et al.||2005] l and 
GRB 0707 14B (Ba rbier et al.T20Q7] i are formally long GRBs, 
but they have spectral features hinting at an underlying coales- 
cence progenitor. GRB 061210 is another long-duration burst, 
but it exhibits th e typical short spikes of a short GRB ([Can] 



nizzo et aT}2 006 ). This gives a list of 33 interesting GRBs with 
which to search for an association with gravitational waves 
from compact binary coalescence. 

Around the trigger time of each interesting GRB, we re- 
quired 2190 s of multiply-coincident data. The detectors oper- 
ated with individual duty cycles of 67% to 81% over the span 
of the S5 and VSR1 runs. Where more than two detectors had 
sufficient data, we selected the most sensitive pair based on the 
average inspiral range, because including a third, less sensi- 
tive detector does not enhance the sensitivity greatly. The one 
exception was GRB 070923, described below. In descending 
order of sensitivity, the detectors are HI, LI, H2, and VI. This 
procedure yielded 11 GRBs searched for in Hl-Ll coincident 
data, 9 GRBs in H1-H2, and 1 in H2-L1. 

In addition to these 21, we analyze GRB 070923 because of 
its sky location relative to the detectors' antenna patterns. The 
antenna pattern changes with the location of a source relative to 

a detector and can be expressed by the response y + , 

i n which F^. and F x denote the antenna-pattern functions 
Allen et al.|2005| l. A value of 1 corresponds to an optimal 
ocation of the putative gravitational-wave source relative to 
the observatory, while a value of corresponds to a source 
location that will not induce any strain in the detector. For 
this particular GRB, the optimal antenna response for Virgo 
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is around 0.7, while those for the two LIGO sites are about 
half that (see Table [2|, yielding a comparable sensitivity in 
the direction of GRB070923 for all three of them. Data from 
HI, LI, and VI were analyzed, making this the only GRB 
involving triple coincidences. 

Table [I] lists all 22 target GRBs after applying the selec- 
tion criteria described in this section. Plausible redshifts have 
been published for only three of these GRBs, placing them 
well outside of our detectors' range, but short GRB redshift 
determinations are in general sufficiently tentative to warrant 
searching for all of these GRBs. 

GRB 070201 is also worth special mention. It was already 
analyzed in a high-priority search because of the striking 
spatial coincidence of this GRB with M31, a galaxy only 
^780 kpc from Earth. No gravitational-wave signal was found 
and a coalescence scenari o could be ruled out with >99% 
confidence at that distance ( [Abbott et al.|2008a| ), lendi ng ad- 
ditional support for a soft gamma repeater hypothesis ( |Ofek| 
et al.|200 8 ). However, because of improvements in the analysis 
pipeline, we reanal yzed t his GRB and report the results in this 
paper. See Section 3.1.1 for details. 



2.3. Candidate generation 

We generated candidates using the standard, untriggered 
comp act binary coa l escen ce search pipeline described in de- 
tail in Abbott et al. (2009e I. The core of the inspiral search 
involves correlating the measured data against the theoreti- 
cal waveforms expected from compact binary coalesce nce, a 
technique called matched filtering (Helmstrom 1968). The 
gravitational waves from the inspiral phase, when the binary 
orbit decays under gravitational-wave emission prior to merger, 
are accurately modeled by post-Newtonian approximants in 
the ban d of the detector 's sensitivity for a wide range of binary 
masses ( Blanc het 2006). The expected gravitational-wave sig- 
nal, as measured by LIGO and Virgo, depends on the masses 
(mNS i w comp ) and spins (s*ns > s com p) of the neutron star and 
its companion (either NS or BH), as well as the spatial loca- 
tion (a, 5), inclination angle i, and polarization angle ip of the 
orbital axis. In general, the power of matched filtering depends 
most sensitively on accurately tracking the phase evolution 
of the signal. The phasing of compact binary inspiral signals 
depends on the masses and spins, the time of merger, and an 
overall phase. 

We adopted a discrete bank of template waveforms that 
span a two-dimensional parameter space (one for each com- 
ponent mass) such that the maximum loss in signal to noise 
ratio (SNR) for a binary with negligible spins would be 3% 
(Cokelaer 2007 ). While the spin is ignored in the template 
waveforms, we verify that the search can still detect binaries 
with most physically reasonable spin orientations and magni- 
tudes with only moderate loss in sensitivity. For simplicity, the 
template bank is symmetric in component masses, spanning 
the range [2, 40) M in total mass. The number of template 
waveforms required to achieve this coverage depends on the 
detector noise spectrum; for the data analyzed in this paper the 
number of templates was around 7000 for each detector. 

We filtered the data from each of the detectors through each 
template in the bank. If the matched filter SNR exceeds a 
threshold, the template masses and the time of the maximum 
SNR are recorded. For a given template, threshold crossings 
are clustered in time usin g a sliding wind ow equal to the du- 
ration of the template (Allen et al. 2005|. For each trigger 
identified in this way, the coalescence phase and the effective 
distance — the distance at which an optimally oriented and 



optimally located binary, with masses corresponding to those 
of the template, would give the observed SNR — are also com- 
puted. Triggers identified in each detector are further required 
to be coincident with their time and mass parameters with a 
trigger from at least one other detector, taki ng into account 
the co rrelations between those parameters (Robinson et al. 
2008}. This significantly reduces the number of background 
triggers that arise from matched filtering in each detector inde- 
pendently. 

The SNR threshold for the matched filtering step was chosen 
differently depending on which detectors' data are available 
for a given GRB. If data from HI and LI were analyzed, the 
threshold for each detector was set to 4.25, reflecting their 
comparable sensitivity. If data from HI and H2 were analyzed, 
the threshold of the latter detector — the less sensitive of the 
two — was set to 3.5 to gain maximum network sensitivity, 
while the threshold of the more sensitive detector, HI, was set 
to 5.5 since any signal seen in H2 would be twice as loud in 
HI, with some uncertainty. In the single case of analyzing only 
H2-L1 data (GRB 070707) the threshold was 4.25 for LI and 
3.5 for H2, and for the single case of analyzing data with Virgo 
(GRB 070923), the threshold was set to 4.25 for all involved 
detectors (HI, LI, and VI). For comparison, a unifor m SNR 
threshold of 5.5 was used in the untriggered S5 search (Abbott] 
|et al.|2009bp . 

We applied two signal-based tests to reduc e and refine o ur 
trigger sets. First, we computed a \ 2 statistic (Allen 2005 i to 
measure how different a trigger's SNR integrand looks from 
that of a real signal; triggers with large y 2 were discarded. Sec- 
ond, we applied the r 2 veto ( |Rodriguez|2007j ) which discards 
triggers depending on the duration that the % 2 statistic stays 
above a threshold. The SNR and x 2 f rom a single detector 
were combined into an effective SNR ( Ab bott et al.|2 008b). 
The effective SNRs from the analyzed detectors were then 
added in quadrature to form a single quantity p 2 ff which pro- 
vided better separation between signal candidate events and 
background than SNR alone. The list of coincident triggers at 
this stage are then called candidate events. 

2.4. Ranking candidates 

The distribution of effective SNRs from background and 
from signals can vary strongly across the template bank, de- 
pending most strongly on the chirp mass, a combination of the 
two component masses that appears in the leading term of the 
signal amplitude and phase (Thorne 1987 ). For this reason, we 
refine our candidate ranking with a likelihood-ratio statistic, 
which we compute for every candidate in the on-source, off- 
source, and injection trials. In short, we define the likelihood 
ratio L for a candidate c to be the efficiency divided by the 
false-alarm probability. The efficiency here is the probability 
of obtaining a candidate as loud or louder than c (by effec- 
tive SNR) within the same region of template space given a 
signal in the data. The efficiency is a function of the signal 
parameters TO comp and D and is marginalized over all other 
parameters; it is obtained by simply counting across injection 
trials. The false-alarm probability here is the probability of 
obtaining a candidate as loud or louder than c in the same re- 
gion of chirp mass from noise alone; it is obtained by counting 
across off-source trials. 

At the end of the search (i.e., Table [2] and Figure [T]), we 
report a different false-alarm probability. It is the fraction of 
off-source likelihood ratios larger than the largest on-source 
likelihood ratio. 

There is another noteworthy difference with respect to un- 
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Table 2 

Summary of the results for the search for gravitational waves from each GRB. 

The Antenna Resp onse column contains the response for each detector as 
explained in Section [2~2) an ellipsis (...) denotes that a detector's data were 

not used. F.A.R is the false-alarm probability of the most significant 
on-source candidat e for a GRB as measured against its off-source trials, as 
explained in Section [2~4] On-source trials with no candidates above threshold 
are assigned a F.A.R of 1 . The last two columns show t he low er limits at 90% 
CL on distances, explained in Section |3~2] 

triggered inspiral searches. For background estimation, untrig- 
gered searches use coincidences found between triggers from 
different detectors, to which unphysical time-shifts greater 
than the light-travel time between detector sites have been 
applied. Unfortunately, HI and H2, being co-located, share a 
common environmental noise that is absent from the time-shift 
background measurement. Being unable to estimate the signif- 
icance of H1-H2 candidates reliably, the untriggered search 
examines them with significantly greater reservation and does 
not consider them at all in upper limit statements on rates. 
The present search performs its background estimation with 
unshifted coincidences under the assumption that any gravita- 
tional wave signal will appear only in the on-source trial. Thus, 
we regain the unconditional use of H1-H2 candidates. 

3. RESULTS 

3.1. Individual GRB results 

We found no evidence for a gravitational-wave signal in 
coincidence with any GRB in our sample. We ran the search 
as described in the previous section and found that the loud- 
est observed candidates in each GRB's on-source segment is 
consistent with the expectation from its off-source trials. The 
results are summarized in Table |2] with brief highlights in the 
following subsections. A graphical comparison of on-source 
to off-source false-alarm probability is shown in Figure [T] 

3.1.1. GRB 070201 

The reanalysis of GRB 070201 yielded candidates in the 
on-source segment, despite havin g no coincident c andidate 
at all in the previous analysis (Ab bott et al.|2008a i. This is 
consistent because the threshold for H2 has been lowered from 
4.0 to 3.5 and the coincident trigger found in this reanalysis 
happened to lie very close to the larger threshold in the previous 
search. The reanalysis yields a false-alarm probability of 6.8%, 
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Figure 1. Cumulative false-alarm probabilities for the most sign ificant candi- 
date in each on- and off-source trial, as described in Section |2l4] 

the smallest in the set of analyzed GRBs l . This value is 
completely within our expectations when we consider that we 
examined 22 GRBs. 

3.1.2. GRB 070923 

GRB 070923 was the GRB for which HI, LI, and VI had 
comparable sensitivity and we accepted triggers from all three 
detectors. There were no triply-coincident candidates in the 
on-source trial, but there were surviving doubly-coincident 
candidates, the loudest of which had a false-alarm probability 
of 74.5%. 

3.2. Distance exclusions 

With our null observations and a large number of simula- 
tions, we can constrain the distance to each GRB assuming it 
was caused by a compact binary coalescence with a neutron 
star (with a mass in the range [1, 3) M ) and a companion 



of mass Wcomp- For a given m comp range, we used the ap- 
proach of Feld man & Cousins] ( |1998| l to compute regions in 
distance where gravitational-wave events would, with a given 
confidence, have produced results inconsistent with our ob- 
servations. Figure [2] shows the lower Feldman-Cousins dis- 
tances for the 22 analyzed GRBs at 90% confidence for two 
illustrative choices for the companion mass range. The val- 
ues are also listed in Table [2j Because the companion mass 
range has been divided into equally spaced bins, we report 
on a 'NS-NS' system in which the companion mass is in the 
range [1, 4) M Q and a 'NS-BH' system in which the BH has a 
mass in the range [7, 10) M Q . The median exclusion distance 
for a NS-BH system is 6.7 Mpc and for a NS-NS system is 
3.3 Mpc. These distances were derived assuming no beaming 

1 In public presentations of preliminary results, GRB 061006 was erro- 
neously highlighted as having the loudest c andidate due to a 22 .8 s offset 
in the GR B time. Swift's initi al GCN alert (Schady et al. 2006a) was later 
corrected I Schady et al. 2006b ), but we initially overlooked this correction. 
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Figure 2. Lower limits on distances at 90% CL to putative NS-NS and 
NS-BH progenitor systems, as listed in Tableland explained in Section |3~2] 

(uniform prior on cos i). NS-BH distances are typically higher 
than NS-NS because more massive systems radiate more total 
gravitational-wave energy. The excluded distance depends on 
various parameters: the location of the GRB on the sky, the 
detectors used for the GRB, the noise floor of the data itself, 
and the likelihood ratio of the loudest on-source candidate 
event for the GRB. 

We drew the simulations from a distribution in which 
our marginalized parameters roughly reflect our pri- 
ors on these astrophysical compact binary systems. 
In our models, a signal is completely specified by 
(m NS ,m comp ,SN S ,s comp ,i,4>,t ,D,a,5). Of these, 
we wish to constrain m com p and D, marginalizing over 
everything else. We drew the NS mass tons uniformly 
from [1,3) Mq\ the magnitudes of t he NS spins Jsns I were 
half and half uniform in [0, 0.75) ( jCook et al J 1994b ; the 
magnitude of the compan ion's spin |s"c 0mD | were hal ft) and 
half uniform in [0, 0.98) ( |Mandel & O'Shaughnessy|2009| ); 
the orientations of the spins were uniform in solid angle; 
the inclination i of the normal to the binary's orbital plane 
relative to our line of sight was conservatively chosen to be 
uniform in cos t instead of making an assumption about the 
GRB beaming angle; the polarization angle <fi was uniform 
in [0, 2ir); the coalescence time to was uniform over the 
off-source region; the declination S was set to that of the GRB; 
the right ascension a was also set to that of the GRB, but 
was adjusted based on to to keep each simulation at the same 
location relative to the detector as the GRB. 

A number of systematic uncertainties enter into this analy- 
sis, but amplitude calibration error and Monte-Carlo counting 
statistics from the injection trials have the largest effects. We 
multiplied exclusion distances by 1.28 x (1 + <5 ca i), where S ca \ 
is the fractional uncertainty (10% for HI and H2; 13% for LI; 
6% for VI ( jMarion et aLl[2008] |). The factor of 1.28 corre- 
sponds to a 90% pessimistic fluctuation, assuming Gaussianity. 
To take the counting statistics into account, we stretched the 



Feldman-Cousins confidence belts to cover the probability 
CL + 1.28^/CL(1 - CL)/n, where CL is the desired confi- 
dence limit and n is the number of simulations contained in 
the (m comp , D) bin for which we are constructing the belt. 

3.3. Population statement 

In addition to the individual detection searches above, we 
would like to assess the presence of gravitational-wave signals 
that are too weak to stand out above background separately, 
but that are significant when the entire population of analyzed 
GRBs is taken together. We compare the cumulative distribu- 
tion of the false alarm probabilities of the on-source sample 
with the off-source sample. The on-source sample consist of 
the results of all 22 individual searches, including those for 
GRBs with known redshifts, and the off-source sample con- 
sists of 6801 results from the off-source trials. This number is 
lower than 22 x 324 because for some GRBs, some off-source 
trials were discarded due to known data quality issues. 

These two distributions are compared in Figure [T] To deter- 
mine if they are consistent with being drawn from the same 
parent distribution, we emplo y the non-parametric Wilcoxon- 
Mann- Whitney U statistic ( |Mann & Whitney|1947) , which is 
a measure of how different two populations are. Applying the 
U test, we find that the two distributions are consistent with 
each other; if the on-source and off-source significances were 
drawn from the same distribution, they would yield a U statis- 
tic greater than what we observed 53% of the time. Therefore 
we find no evidence for an excess of weak gravitational-wave 
signals associated with GRBs. 

4. DISCUSSION 

We searched data taken with the three LIGO detectors and 
the Virgo detector for gravitational-wave signatures of compact 
binary coalescences associated with 22 GRBs but found none. 
We were sensitive to systems with total masses 2 Mq < to < 
40 M Q . We also searched for a population of signals too weak 
to be individually detected, but again found no evidence. While 
there are few redshift determinations for short GR Bs, it appears 
that the distribution is peaked around (z) ~ 0.25 (Naka r|2007} , 
far outside initial detector sensitivity, so it is not surprising that 
the S5/VSR1 run yielded no detections associated with short 
GRBs. 
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